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Abstract 

In this paper we study the connection between four models describing dislocation dynamics: a generalized 2D Frenkel- 
Kontorova model at the atomic level, the Peierls-Nabarro model, the discrete dislocation dynamics and a macroscopic 
model with dislocation densities. We show how each model can be deduced from the previous one at a smaller scale. 
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1 Introduction 

In this paper, we focus on the modelling of dislocation dynamics. We refer the reader to the book of Hirth and 
Lothe [20] for a detailed introduction to dislocations. Our study ranges from atomic models to macroscopic 
crystal plasticity. At each scale, dislocations can be described by a suitable model. Our goal is to explain how 
we can deduce a model at a larger scale, from the model at a smaller scale. 

Even if our derivation will be done on some simplified models (essentially 2D and ID models), we hope 
that our contribution will shed light, even on some well-known models. More precisely, we will consider the 
following four models, from the smaller to the larger scale: 

1. Generalized Frenkel-Kontorova model (FK) 

2. Peierls-Nabarro model (PN) 

3. Dynamics of discrete dislocations (DDD) 

4. Dislocation density model (DD) 

Schematically, the four models are related as shown below (see also Figure 6.1 for a more detailed diagram): 



(FK)e, 



(PN) 



e2,e3>0 



£2^0 



(DDD) 



e3>0 



£3^0 



(DD) 



(1.1) 



The rest of the paper is composed of four sections. Each section presents one model, and explains how this 
model can be deduced from the previous model at a smaller scale. 



2 Generalised Frenkel-Kontorova model 

2.1 Geometrical description 

We call (61,62,63) a direct orthonormal basis of the threedimensional space. We consider a perfect crystal 
Z'^ where each position with integer coordinates is occupied by one atom. We want to describe dislocations, 
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which are certain "line defects" in the crystal. To simplify the presentation, we will assume that the material 
is invariant by integer translations in the direction 63. Because of this assumption, we can simply consider the 
cross section of the crystal in the plane (ei, 62) where each atom is now assumed to have a position / G in 
the perfect crystal. We also assume that each atom / can have a displacement [// G R in the direction ei, such 
that the effective position of the atom J is / + Ujei. 

On Figure [1] below is represented a view of the perfect crystal. On Figure [2] we can see a schematic view 
of a edge dislocation in the crystal. On this picture, the upper part {I2 > 0} of the crystal has been expanded 
to the right of a vector ^ei, while the lower part {I2 < —1} of the crystal has been contracted to the left of 
a vector ~\ei- The net difference between these two vectors is ei and is called the Burgers vector of this 
dislocation. 



I2=0 
l2=-l 



^1 

Figure 1: Perfect crystal 



l2=0 
l2=-l 



Figure 2: Schematic view of a edge dislocation in 
the crystal 



In order to describe a edge dislocation in our formalism, let us make a few assumptions. We will assume that 
the dislocation defects are essentially described by the mismatch between the two planes I2 = Q and I2 = — 1, 
like on Figure [H For this reason, and also in order to simplify the analysis, we assume that the displacement 
of the crystal satisfies the following antisymmetry property 

Uiiu-h) = -Uih.i,-i) for all /=(/i,/2)eZ2. (2.1) 

Let us also define the discrete gradient 

Remark that defects in the crystal can be seen as 



Ul+e, - Ul 
Ul+e, - Ul 

regions where the discrete gradient is not small. 



Formalism for a edge dislocation with Burgers vector ei 

In our formalism, a edge dislocation like the one of Figure [21 can be represented by a displacement Ui satisfying 

C^(/i,o) ^ -1) ^0 as /i ^ -00 

U{i,,o) = -C^(/i -1) ^ ^ as /i ^ +00. 

Because we assume that the dislocation core lies in the two planes I2 — and I2 = —1, it is reasonable to 
assume that all the components of the discrete gradient are small, except components C//+e2 ^ Uj for / = (/i, I2) 
with I2 = —1. More precisely, we assume that there exists a small 5 > such that 



\UI+e^ -Ui\<S for aU / = {h,^) G 
\Ui+e2 -Ui\<S for aU / = ih,^) G with I2 ^ -1 



2 „ruu -I 1 (2-2) 
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Moreover, if there is no apphed stress on the crystal, then it is reasonable to assume that 

dist ((V''C/)/,Z2) ^ as |/|^+oo. 

2.2 The energy and the dynamics 

We assume that the energy of a configuration U = {Ui)i(zz of the crystal can be formally written as 

where 14^ : M ^JR is a potential describing nearest neighbors interactions satisfying 
Assumption (Al) 

fi) (Regularity) W e C^{R) 

ii) (Periodicity) W{a + 1) = W{a) for aU a £ R 

iii) (Minimum on Z) W{Z) = < W{a) for aU a e M\Z 

iv) (Local harmonicity of W) W{a) — for all \a\ < S 

where (5 > is introduced in (|2.2p . Remark that the periodicity of the potential W reflects the periodicity 
of the crystal, while the mimimum property of W is consistent with the fact that the perfect crystal is 
assumed to minimize its energy. Assumption iv) will be used for later simplification. 

Then we assume that we are in a regime where the crystal reaches very quickly the equilibrium in the 
regions where there is no defects, i.e. satisfies 

= -\7ujE{U) for all / = (h^h) e Z^ with I2 ^ 0, -1 (2.3) 

while we have the following fully overdamped dynamics in the two planes where the dislocation lives (describing 
the average friction of the lattice on the effective dissipative motion of the dislocations): 

^Ui = -VurEiU) for aU / iluh) e 1? with h = 0, -1. (2.4) 
dt 

Let us mention that we do not have a fundamental justification of this dynamics, but we think that one of the 
main justification of this model is that other known models at larger scales can be deduced from this particular 
model. For physical justifications of the dissipative effects in the motion of dislocations, see [21 [20]. See also 
[22| 1231 124] for a fundamental justification of the overdamped dynamics based on explicit computations in a 
ID Hamiltonian model. 

Taking into account the local harmonic assumption (Al) iv), applied where the components of the discrete 
gradient are small (see (|2.2I) ). joint to the antisymmetry property defined in (|2.ip . we can rewrite system 
([Ml)- ([23) as follows for all t > 0: 

'0= E (Uj-Ui) foraU I ^ {h, h) with I2 > I 

JGZ^, |J-/| = 1 

— Ui = -W'(2Ui)+ V (Uj-Ui) foraU I = {h, h) ^ '^^ with /s = 0. 

(it ^ 

|J~/| = 1, J2>0 

We call this model a generalised Frenkel-Kontorova model. Even if this system of equations is not standard, 
it is nevertheless possible to define a unique solution under suitable assumptions, in the framework of viscosity 
solutions (see [IS])- We refer the reader to the book of Braun, Kivshar, [5] for a detailed presentation of the 
classical FK model. For homogenization results of FK models, we refer the reader to [15]. For the description 
of vacancy defects at equilibruim, see [19]. See also [18], where the authors study the problem involving a 
dislocation inside the interphase between two identical lattices. Their model corresponds to our model (|2.5p at 
the equilibrium with the potential is a cosinus function. For other 2D FK models, see [8] [9]. 
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Remark 2.1 It is important to remark that we used condition i2.2\) only to derive the model. We do not know 
and we do no claim that there exists solutions of system \2.5]) satisfying condition \2.2\) . From now on, we 
only consider solutions of system \2.5]) without requiring further assumptions on the solutions. 

Remark 2.2 When we freeze the components Ui — for I2 > 1, and change the evolution equation forgetting 
the index J with J2 — I, this leads to the following classical fully overdamped Frenkel-Kontorova model satisfied 
by Vi := C/(i,o) 

2.3 The asymptotic stress created by a single dislocation 

In this subsection, wc will compute the asymptotic stress created by a single dislocation. To this end, we first 
compute the effective Hook's law of the lattice. 
Computation of the Hook's law 
Let us consider an affine displacement 

Ui ^ a ■ I + C with a = (ai , a2) e 

where C G M is a constant. Then the energy by unit cell is 

£ = WiUi+e, - Ui)) + WiUi+e, - Ui)) = lial+ al) 

for \a\ < 6. Reminding the fact that U is the displacement in the ei direction, we get that the strain e (i.e. 
the symmetric part of the gradient of the displacement) is given by 

''')-lm^e,+e,^yU)^(^\^ 

621 622/2 V 02/2 



Therefore 



^(6) = ^e2i+2e?2. 



d£ 

Recalling that the stress is given by cr° = — , we get the Hook's law: 

oe 

_ f 611 2ei2 

^ ^ \ 2621 

Computation of the stress created by a single dislocation 

Remark that when there is no dislocations, the energy associated to a continuous displacement U{X) for 
X ^{Xi,X2) is formally 

E = [ J|V[/p. 

Therefore the Euler-Lagrange equation (which is the corresponding equation of elasticity for this model) is 

AU = 0. (2.6) 

Let us now consider the following function 

Uo{X) = ^ arctan + ^sgn (X^) 

where sgn is the sign function. This function satisfies 

f UoiXi,X2) = -UoiXi,-X2) 

Uo{Xi,0+) = -UoiXi,0-) ^0 as Xi ^ -00 
i7o(Xi,0+) = -C/o(Xi,0-)^i as 
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Moreover we can easily check that 

div (yUo- H{Xi)So{X2)e2) ^0 in ^'(M^) 

where H is the Heavyside function and is the Dirac mass. This equation is the analogue of equation (|2.6p 
when there is a dislocation. This shows that in a continuum mechanics framework associated to the particular 
lattice that we consider, the function Uq is the displacement corresponding to a dislocation with Burgers vector 
ei. In particular, the stress created by this dislocation is then given by 



1 



xl + xl xl + xj 



and then 



aO2(Xi,0) = -^ (2.7) 

which is the asymptotic shear stress at the point {Xi , 0) G created by a single dislocation positioned at the 
origin, and with Burgers vector ei. 

2.4 Rescaling of the generalised FK model 

In this subsection, in order to simplify the notation we denote by e := £i > the small parameter in the first 
passage of the scheme (|l.ip . We are interested in the case of asymptotically small potential W for which we 
expect an asymptotically large dislocation core. This means that in this limit, we expect to be able to describe 
the discrete displacement Ui by a continuous function. 

More precisely, we first define the rescaled integer coordinates: 

Vf" = (eZ) X e (N\ {0}) , dVL" = (eZ) x {0} . 

Then we write the potential as 



W = 



and define the rescaled function 



u''{X,t)^2Ux{-] for X={Xi,X2)en\ te[0,+oo). 



Remark that the factor 2 in the definition of permits to interprete as the jump of the displacement in 
the direction ei, when we pass from hyperplane X2 = — e to the hyperplane X2 = 0. 

We can easily check that solves the following system of equations (with the particular value a = 0) 

= 4 (u^X + eJ,t)-u%X,t)) for all {X,t) £ x {0,+oo) 

^ jez2, |j|=i 

ul{X,t) ^ 662(7 (£62X1) ~ (Wyiu^X^t)) + F[u%X,t) foraU (X, t) e 51)^ x (0, +00) (2.8) 
with r[u'']{X,t) = - {u-{X + eJ,t) ~u''{X,t)). 

^ |J| = 1, J2>0 

Here £2 > is a small parameter, and the scalar function a has been introduced to take into account the 
possible external applied shear stress on the material. We will also assume that the initial data satisfies 

u^X,0) ^uoiX) foraU X e dn" (2.9) 
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where uq is a given function independent on s and smooth enough. 

In order to identify a hmit model as e goes to zero, we also make the following assumption 



\\W'-W\\ci 











(2.10) 



for some new potential W satisfying the following assumption: 
Assumption (Al) 

The potential W satisfies (^1) i), ii), iii) 
iv) (Non degenerate minima): a := W"{Q) > 0. 



In (|2.10p . we use the norm, because this is the first derivative of the potential that appears in the 
equations. Remark that condition (j2.10p can be fulfilled, if we assume for instance that W satisfies assumption 
(Al) with S ^ << e. We also make the following assumption on the stress: 
Assumption (A2) 

There exists a constant C > such that 

W\ + Wx\ + Wxxl < C for all x e M. 



3 The Peierls-Nabarro model 
3.1 Description of the PN model 

In this section we introduce the Peierls-Nabarro model, which is a phase field model (see 
of this model). In this model, phase transitions describe the dislocation cores. We set 



for a presentation 



{X ^ {Xi,X2) X2>o}. 

A function u°(X, t) is said to be a solution of the PN model, if it satisfies the following system 



u'l ^2s2a{s2Xi)^W'iu°) + 



dvP_ 
dX^ 



on Q X (0, +oo) 
on X (0, -l-cxo) 



(3.1) 



with initial data 



i°(X, 0) = uo{X) for all X e dfl. 



(3.2) 

The stationary version of this model has been originally introduced in order to propose a method to compute 
at the equilibrium a finite stress created by a dislocation. In this model, is the phase transition. For instance, 
for a edge dislocation with Burgers vector ei as presented in Section [21 is a transition between the value 
on the left to the value 1 on the right (see Figure [3]) . In the special case = = cr and for sinusoidal 
potentials W, the stationary solution u° is known explicitely (see for instance [7 ), which makes the PN model 
very attractive. Let us mention that a physical and numerical study of the evolution problem (|3.ip has been 
treated in 1^. 



1 -- 







u 







Figure 3: Phase transition for a edge dislocation with Burgers vector ei for X2 — 
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Remark 3.1 Remark that when we consider system (K7\l-(^^ in the framework of viscosity solutions, the 
evolution equation on the boundary dfl appears to be a boundary condition of the system. For this reason, as it 
is usual for viscosity solutions (see for instance this boundary condition has to be understood technically 

in the sense that on dfl the function u'^ solves pointwisely either — or the evolution equation. 



3.2 Convergence of the generalised FK model to the PN model 

We have the following result 

Theorem 3.2 (Formal convergence of FK to PN, e = ei ^ 0) 

Let £ > 0. For the initial data uq G W'^'°°{fl) which is assumed harmonic on fl, and under assumption (Al) 
on eW^ , and (A2) on a, there exists a unique viscosity solution of system h2.8\) - l27^) . Moreover assuming 
112.1 0\) with the potential W satisfying assumption (Al), then, as e goes to zero, the solution formally 
converges to a viscosity solution of system Iff. Jp -i TO)) . 

The proof of Theorem 13.21 is done in full details in [T3] . 
Sketch of the proof of Theorem 13.21 

One way to guess the limit model (|3.ip - (|3.2p is to pass to the limit formally in system (|2.8p - (|2.9p assuming that 
the solution u'^ (and its derivatives) converges to a limit . The convergence in the system is then obtained 
using a simple Taylor expansion. The existence of a solution to system (I2.8|) - (|2.9p is technically delicate 
and is based on the proof of a suitable comparison principle for this system. 



3.3 Reformulation of the PN model 



In this subsection and in Subsection 13.41 in order to simplify the notation we denote by e := £2 > the small 
parameter in the second passage of the scheme (jl.ip . We recall that it is well known that for bounded smooth 
functions defined on J7 which are harmonic on SI, we can write 



dX~2 



(Xi,0) = L(u"(-,0))(Xi) for all {Xi,0) e dQ 



where for a general bounded smooth function w, the linear operator L is given by the Levy-Khintchine formula 
(see Theorem 1 in PTj): 



iLw)ix) 



dz 



{■w[x + z) - w{x) - 2;w'(a;)l{|^|<i}) . 



(3.3) 



Then for smooth solutions system (j3.1l) can be rewritten for V{X\,t) = vf^{X\,X2,t)\XQ,=a with x = 
Xi e R as 

Vt = 2ea{ex) - W' (V) + LV on M. (3.4) 
We also recall (see [7[) that there exists a unique function solution of 

' L<()-W'{^) on M 
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(3.5) 



-(X3) = 0, 0(0) = -, 0(+^) = l. 



0' > and 

The function is called the layer solution and a translation of is pictured on Figure [31 



3.4 Rescaling of the PN model 

We now consider the following rescaling 
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Lv' - -W'iv") + 2a{x) } on R (3.6) 



Then system p.4p can be rewritten as 

1 



with initial condition 

v%x,0) ^v^q{x) for xeR. (3.7) 

Again, a good notion of solution for system (|3.6p - (|3.7p is the notion of viscosity solution for non local 
equations (see for instance [5]). 

Here we will choose carefully the initial condition Vq as follows 
Assumption (A3) 

Ju-y \ ••■ -^jY 

^ / 

£ , ^ ■^-^ , X — Xi 



where we recall that a = W" (0) > and (j) is defined in 



4 Dynamics of discrete dislocations 
4.1 Description of the DDD model 

In this section we assume that the phase transition reduces to a sharp interface where the transition is localized 
at the position x = x\ ^M.. For a dislocation associated to a Burgers vector ei, the sharp interface is associated 
to a non-decreasing step function like the one of Figure |4l 



1 + 





Figure 4: Sharp interface describing a discrete dislocation at a; = a; 5 

More generally, we can consider the case of dislocations (or particules) of positions {xi{t))i^i^,,,^N solving 
the following system 

'^ = ^^\a{x,)+Y,y'i^^-x,)] on (0,+oo) for i = l,...,A (4.1) 
with the two-body interaction potential 

V{x)^~^\n\x\ 

with initial data 

x^[0)^x° for i = l,...,A. (4.2) 

Here the constant 7 > is the inverse of the damping factor. It is related to the layer solution defined in 
(|3.5p and is given by 

-1 

' 



7 = 2 (0') 
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The function a is the apphed shear stress and V {x — Xj) is the (singular) shear stress created at the point 
X by the dislocation xj. This corresponds exactly to the shear stress already computed in (|2.7p . The total 
stress <j{xi) + V'{ called the resolved Peach-Koehler force acting on the dislocation Xi. 

4.2 Convergence of the PN model to the DDD model 

We have 

Theorem 4.1 (Convergence of PN to DDD, £ = £2—^0) 

Let £ > 0. Under assumptions (A1)-(A2)-(A3), there exists a unique viscosity solution v'^ of if5'.6]j - |X7[ ). 
Moreover there exists a unique solution of Jp - JX^ , and we define 

v\x,t)^ J2 H{x~x,it)) 

i=l,...,N 

where H is the Heavyside function. Then as e goes to zero, the function converges to in the following 
sense 

limsup v%x',t') < {v°y{x,t) 

{x' ,t')~>{x,t), e^O 

and 

liminf v^ix^t') > {v°)Jx,t). 

{x',t')^{x,t), e^O 

The proof of this result is done in full details in [T7] . 
Remark 4.2 We recall that the semi- continuous envelopes of a function v are defined as 
v*{x,t)= limsup v(x',t') and v<,{x,t) — liminf v{x',t'). 

{x',t')-^(x,t) (x',t')^{x,t) 

Sketch of the proof of convergence 

The existence of a solution for all time of the ODE system ()4.ip - ()4.2|) comes from the fact that V{x) is a convex 
potential outside the origin. This property allows to show that the minimal distance between particles 

d{t) = inf \x,{t) ~ x.j{t)\ 

satisfies 

d{t) > d(0)e-^'^* (4.3) 
which prevents the meeting of the particles at any finite time. 

Then the main idea to prove the convergence is to approximate the solution t;^ by the following ansatz 

\ 1 , / s dxj , , 

w^(x,i) = --2a(x) + 2]|0(^^^j -e2;,(i)V'(^^^j| with x,{t) ^ -^{t) 
where a = W"{0) and the corrector ip solves the following equation 

Lij-W"i^)^^^' + v{W"i<P)~W"m with rj = —L_Jy)^ 
The stress created in a; by a dislocation positioned at the origin, comes from the following property 

(j)(x) — H(x) ~ as \x\ +0O. 

aiTx 

The rest of the proof of convergence of is done by construction of sub and super solutions based on the 
ansatz v^. 
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4.3 Rescaling of the DDD model 

In this subsection, in order to simplify the notation we denote by e := £3 > the small parameter in the third 
passage of the scheme (jl.ip . We consider a given initial data wo which satisfies 
Assumption (A4) 

Wo £ VF2'°°(R), 
w'q > 0, wq{—co) = 0. 

We also introduce the integer A^^ and the position of the dislocations Xi < ... < such that 

i=l,...,N^ 



where [-J denotes the floor function. We also assume that the stress a is periodic. Precisely, we make the 
following assumption: 
Assumption (A2') 

<T e C^(R) and a{x + 1) = (j{x) for all a; e R. 

This assumption allows to study the collective behaviour of dislocations in a landscape with periodic obstacles, 
and to get the effective macroscopic model by a periodic homogenization approach. 

Then we consider the solution {xi{t))i^i jy^ of the system (|4.ip - (l4.2p with N — and define the function 



v°{x,t)= H{x-x,{t)) 



i=l,...,Nc 

and the rescaling 

w'{x,t)^ev'' (^,-^ . (4.4) 

5 Dislocation density model 

5.1 Description of the DD model 

We first introduce a function g : (0, +00) x M ^ M which satisfies 
Assumption (A5) 

" g e C°((0,+cx)) X M), 

I I— > g[p, I) is nondecreasing. 
Then we consider a function w^{x, t) which is a solution of 

w^t = g{wl, LuP) on R x (0, +00) (5.1) 

where the operator L is defined in (|3.3p . and with initial data 

w°{x, 0) ^ wo{x) for all a; G R. (5.2) 

Here the function vP is such that its derivative represents the macroscopic dislocation density. Moreover 

can be seen as the plastic strain localized in plane X2 — Q and i^LuP can be identified to the stress created by 

the dislocation density w° . Equation (|5.ip can be interpreted as the plastic flow rule in a model for macroscopic 
crystal plasticity. Indeed, from a mechanical point of view, we have the following table (see also [21]) of 
equivalence between our homogenized model and a classical model in mechanics for elasto-visco-plasticity of 
crystals (see [IB]). 
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Crystal elasto-visco-plasticity 


DD model 


resolved plastic strain 




^i'"(a;l)'5o(a;2) 


Hisloea tiorm Hermities 


a = (ei (g) e2h'{xi)So{x2) 


a = (ei (g e2)w°ixi)So{x2) 


exterior applied stress 


Y^ext 




microscopic resolved 
shear stress 




"-J 


f a 

(0,1) 


resolved exterior 
applied stress 




J 


f a 

(0,1) 


displacement 


V = vid 




strain 


e := e{v) - e^'jSo{x2) 

with e{v) -.^ i {Vv + *Vw) 
1 2 

and e" := - (62 g) ei + ei (g) 62) 

Z 




total elastic energy 


E := [ -{A : e) : e + E'^"'* : e 




macroscopic stress 


E := A : e + E'^^* 




resolved macroscopic 
shear stress 


r := E : e" 




r a 

(0,1) 


visco-plastic law 






2ct ) 

) / 


energy decay 


|i^= /'-r/(r)<0 


^E= f-rg{wl2T)<0 
at Jk 



Remark that when we choose the microscopic stress a in Assumption (A2') so that / cr = 0, we expect 

^(0,1) 

that g satisfies 

g(p,0) = 0. (5.3) 

This equality reflects the pinning of dislocations (see [HI Th 2.6]). In the model presented in the previous 
table, the plastic strain velocity is prescribed by the function g (whose typical profile can be seen on Figure 
[5]) which is assumed to satisfy (j5.3p . 




Figure 5: Graph of the map I 1^ g{p,l) showing a threshold effect 

Let us mention that a similar profile for g as a function of I only have been obtained rigorously for a different 
model involving the motion of a phase transition solution of AUen-Cahn type equations (see [10]). 
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5.2 Convergence of the DDD model to the DD model 

We have the foUowing residt 



Theorem 5.1 (Convergence of DDD to DD, £ = £3-^0) 

Let us assume (A2'), (A4-). Then there exists a function g satisfying assumption (A5). Moreover the function 
defined in J^.^p converges to the unique solution uP of ^5.1\ )- ^57^) . locally uniformly on R x [0, +00). 

The proof of this result is done in full details in [T^ . 

Remark that Theorem 15. II is an homogenization result in the periodic setting. In the particular case where 
the periodic stress a is equal to zero, we get 

g{p,l)^lpl. 

The presence of a non-zero 1-periodic stress with zero mean value, creates a threshold phenomenon where for 
a fixed dislocation density p, the quantity g{p, I) can be equal to zero if |^| is small enough (see for instance the 
numerical simulations in '12] which look like Figure [5|) . 

Sketch of the proof of Theorem 15.11 

Step 1 : Formal determination of the function g 

To determine the function g, we can look formally for Xi{t) defined for all i G Z, which are particular solutions 
of the ODE system (|4.ip with a replaced by '^^ such that 

Xi{t) ^ h(vt + i/p) with h{a + 1) = 1 + h(a) for aU aGM 

where such a function h is called a hull function. Both h and the constant v have to be determined. It can be 
shown that v is unique. Then we set 

g{p,l) = -vp 

which is known in physics as the Orowan's law. 
Step 2: Regularization at short distances 

To avoid the singularity of the potential V{x) — ln|a;|, we can first approximate it by the following 

27r 

symmetric and continuous potential 



Vsix) 



V{x) if |x| > ,5 

linear if x e {-S, 6)\ {0} . 



We consider the function ti"(a;, t) = n H{x — Xi{t)) associated to the dynamics (I4.ip where the potential 

V is replaced by Vs- Then it is possible to show that w° satisfies the following equation with I = 

= |«°|7 |a(x) + L+Ms[v'i;t)]{x)^ (5.4) 

where for a general function w{x), we can define the non local operator 

Ms[w\{x)^ I dz Vi {z)E {w{x + z) - w{x)) with E{a) = - + k if fc<a<fc + l, fceZ 

J\z\>5 2 

where is a odd integer part function. This is possible to introduce a suitable good notion of viscosity solution 
for equation (|5.4p (see [Sj). In particular, we can show that if v^{x,Q) = px, then v°{x,t)/t — > gs{pj) as 
t +00. Moreover it is possible to show the following estimate 

|g,(p,0"ff(p,OI<-^- (5.5) 
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Step 3: Sketch of the proof of convergence in the regularized case 

After a rescaling of the solution u° of (|5.4p with / = 0, we see that w'^{x,t) — ev'^{x/e,t/e) solves an equation 

wl = \w'J-f {aix/e) + Ml[w'i;t)]ix)} 

for some rescaled non local operator Af|. More generally, any continuous solution of the previous equation, 
can be formally written as 

t) ~ w°{x, t) + er{x/e) 

where r is a suitable corrector. One fundamental remark is that as e goes to zero, we can asymptotically split 
the non local term 

Ml[w'{.,t)+er{./s)]ix) ^{Lw'>){x) + 5[r, u;°(x)](a;/e) 

into its long range contribution ^Lw^ and a short range contribution S involving the corrector r. Remark that 
this long range contribution ^Lw'^ = 1/2 is related to the introduction of the constant 1/2 into equation ()5.4p 
used in the definition of gs{p, I). Taking into account this asymptotical splitting, it is then possible to show the 
convergence of to the solution of (|5.ip - (|5.2p with g replaced by gs- The proof can be done in the framework 
of viscosity solutions, adapting the Evans' perturbed test function method. 
Step 4: Sketch of the proof of convergence in the singular case 

The singular case can be reached using an approximation argument. On the one hand, estimate (|4.3p insures 
that the dynamics (|4.ip on the time interval (0, T/e) is equivalent to the same dynamics with V replaced by 
Vs for 5 < = d{0)e~'''^^^'^ . On the other hand, estimate (|5.5|) is independent on e. Then choosing 6 = 6^, 
the convergence of the solution w'^ on the time interval (0,r) can then be obtained by an adaptation of the 
arguments in the regularized case. 

6 Conclusion 

We considered a two-dimensional Frenkel-Kontorova model in the fully overdamped regime. From this model, 
we derived by a scaling argument the time-dependent Peierls-Nabarro model. Looking at the sharp interface 
limit of the phase transitions of the Peierls-Nabarro model, we were able to identify a dynamics of particles 
that corresponds to the classical discrete dislocation dynamics, in the particular case of parallel straight edge 
dislocation lines in the same glide plane with the same Burgers vector. Considering the motion of these 
particles in a landscape with periodic obstacles, we were able to identify at large scale an evolution model for 
the dynamics of a density of dislocations. This model is a macroscopic model for crystal visco-elasto-plasticity, 
where we predicted a plastic flow rule. This last model shows in particular a threshold effect where dislocations 
can be pinned in the obstacles, if the effective stress acting on these dislocations is too small. 

In order to present a summary of our approach, we give here a diagram (see Figure 6.1) that shows the 
links between the four models treated in this paper. 

Up to our knowledge, this derivation of classical models from a single microscopic model (the 2D Frenkel- 
Kontorova model), seems new and allows to make clear connections between different modelling of dislocation 
dynamics. 
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2D Frenkel-Kontorova (FK) 



' = Yl (Ux+At) - Ux{t)) 

J£Z2, |J| = 1 



for X = {Xi,X2) with X2eN\{0} 



f^Ux{t)=eie2a{eie2X,)-^{Wy{2Ux{t))+l'[Ux{t)] for ^2 = 

where I^[Ux{t)]= ^ Ux+j{t) - Ux{t) 
Jez2,|j|=i, j2>o 



£1 

I 





u^^{X,t)=2U^(^-^^^u°{X,t), as £1^0 
W^^ ^ W, as' £1^0 



Peierls-Nabarro (PN) 
vt{x,t) = 2s2a{e2x) - W'{v{x,t)) + Lv{x,t) 



£2 

i 





z;^==(x,t)=v(^,^) as £2^0 

v%x,t)= ^ 



Dynamics of Discrete Dislocations (DDD) 



for i = !,•••, TV 



£3 

I 





a{x + 1) = ct{x) 



Dislocations Density (DD) 



Figure 6: Descriptive diagram summarizing the linlcs between the four models 
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